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ABSTRACT 

Recent galaxy formation models successfully reproduce the local luminosity function (LF) of galaxies by in- 
voking mechanisms to suppress star formation in low- and high-mass galaxies. As these models are optimized 
to fit the LF at low redshift, a crucial question is how well they predict the LF at earlier times. Here we compare 
recently measured rest-frame V-band LFs of galaxies at redshifts 2.0 < z < 3.3 to predictions of semianalytic 
models by De Lucia & Blaizot and Bower et al. and hydrodynamic simulations by Dave et al.. The models 
succeed for some luminosity and redshift ranges and fail for others. A notable success is that the Bower et 
al. model provides a good match to the observed LF at z ^ 3. However, all models predict an increase with 
time of the rest-frame V-band luminosity density, whereas the observations show a decrease. The models also 
have difficulty matching the observed rest-frame colors of galaxies. In all models the luminosity density of red 
galaxies increases sharply from z 3 to z ~ 2.2, whereas it is approximately constant in the observations. Con- 
versely, in the models the luminosity density of blue galaxies is approximately constant, whereas it decreases 
in the observations. These discrepancies cannot be entirely remedied by changing the treatment of dust and 
suggest that current models do not yet provide a complete description of galaxy formation and evolution since 
z~3. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: fundamental parameters — 
galaxies: high-redshift — galaxies: luminosity function, mass function 



1. INTRODUCTION 

In the current paradigm of structure formation, dark mat- 
ter (DM) halos build up in a hierarchical fashion through the 
dissipationless mechanism of gravitationally instability. The 
assembly of the stellar content of galaxies is instead gov- 
erned by much more complicated physical processes, often 
dissipative and non-linear, which are generally poorly under- 
stood. To counter this lack of understanding, prescriptions 
are employed in the galaxy formation models. One of the 
fundamental tools for constraining the physical processes en- 
coded in these models is the luminosity function (LF), since 
its shape retains the imprint of galaxy formation and evolution 
processes. 

The faint end of the LF can be matched with a combina- 
tion of supernova feedback and the suppression of gas cool- 
ing in low-mass h alos due to a backg round of photoioniz- 
ing radiation (e.g.. iBenson et al."1l2002l) . Matching the bright 
end of the LF has proven more challenging. Very recent im- 
plementation of active galactic nucleus (AGN) feedback in 
semianalytic models (SAMs) has yielded exceptionally faith- 
ful reproductions of the observed local rest-frame B- a nd K- 
band global LFs (Bower et al. 20061 ICroton et al. Il200 6: see 
also Granato et al. 2004) , including good matches to the lo- 
cal rest-frame B-band LFs of red and blue gala xies (although 
with some discrepancies for faint red galaxies; ICroton et al. I 
|2006). 

The excellent agreement between observations and mod- 
els at z '~ is impressive but is partly due to the fact 
that the model parameters were adjusted to obtain the best 
match to the local universe. A key question is there- 
fore how well t hese models p r edict t he LF at earlier times. 
The SAMs of Croton et al. (2006]), iDe Lucia & Blaizot I 
(12007) . and B ower et al. ( 200a) have been compared 
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to observations at < z < 2 (see iBower et al. I 120061: 
iKitzbichler & White "2007). Although the agreement is gen- 
erally good, Kitzbichler & White ( 2007) infer that the abun- 
dance of galaxies near the knee of the LF at high redshift is 
larger in the SAMs than in the observations (except possibly 
for the brigh test objects), in an apparent reversal of previous 
studies (e.g.. lCimatti et al. Il200 2). 

Recently, the rest-frame optical LF has been accurately 
measured in the redshift range 2.0 < z < 3.3, using a combi- 
nation of the /:-selected MU SYC, GOODS, and FIRES sur- 
veys (iMarchesini et al.l[2007l) . In this Letter we compare the 
observed LF to that predicted by theoretical models in this 
redshift range, in order to test the predictive power of the lat- 
est generation of galaxy formation models. We also compare 
the observed LF to predictions from smoothed particle hy- 
drodynamics (SPH) simulat ions, which have s o far only been 
compared to data at z ~ 6 (iDave et al. 1 120061) . We note that 
these comparisons ar e effectively the rest-frame equ ivalent of 
the test proposed by iKauffmann & Chariot! (Il998l) . We as- 
sume Q,u = 0.3, = 0.7, and i/o = 70 km s"' Mpc"'. All 
magnitudes are in the AB system, while colors are on the Vega 
system. 

2. THE OBSERVED LUMINOSITY FUNCTIONS 

The obs erved rest-frame optical LFs at z > 2 have been 
taken from IMarchesini et all (l2007h . Briefly, they presented 
the galaxy LFs in the rest-frame B-band (at 2.5 < z < 3.5 
and 2 < z < 2.5), V-band (at 2.7 < z < 3.3), and R-hand (at 
2 < z < 2.5), measured from a /T-selected sample constructed 
from the MUltiway elength Survey by Yale-Chile (MUSYC; 
lOuadri et al. 1120071) . the ultradeep Faint InfraRed Extragalac- 
tic Survey (FIRES; i Franx et al. I l2003h . and the Great Ob- 



servatories Origins Deep Survey (GOODS; iGiavalisco et al.l 
l2004l Chandra Deep Field-South). This /iT-selected sam- 
ple, comprising a total of ^-^990 galaxies with Kl°^ < 25 at 
2 < z < 3.5, is unique for its combination of surveyed area 
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(^380 arcmin^) and large range of luminosities. 

In this Letter we limit our comparison between observed 
and predicted LFs to the rest-frame V band, a t the two redshift 
interv als 2.7 < z < 3.3 (directly taken from Marchesini et alj 
2007) and 2 < z < 2.5 (derived in the same way as described in 
Marchesini et al . 2007). The results are qualitatively similar 
for other rest-frame bands. 

3. THE MODEL-PREDICTED LUMINOSITY FUNCTIONS 

The lBower et al. I (12006 ') SAM is implem ented on the Mil - 
lennium DM simulation described in Sp ringel et al. I (l2005h . 
The details of the assumed prescript ions and the spe - 
cific paramete r choices a re described in ICole et al.1 (120001) . 
iBenson et all (12003). and lBower et al. I (l2006l). We have a lso 
used the outputs^ from the SAM of Cr oton et al. I ^06) as 
updated by De Lucia & Blaizot (2007). This model differs 
from the SAM of iBower et al. I (l2006l) in many ways. The 
scheme for building the merger trees is different in detail, 
as are many of the prescriptions adopted to model the bary- 
onic physics, most notably those associated with the growth of 
and the feedback from SMBHs in galaxy nuclei and the cool- 
ing model (see|K auffmann & Haehnelt 2000; Springe l et al. I 
l200UlDe Luciaet al. .2004: De Lucia & Blaizot ■2007i for de- 
tails). Finally, we have compared the observed LFs with 
the predictions from the cosmological SPH simulat ions of 
lOppenheimer & Dave (2006), already used in Finlator et al. I 
(120078) to constrain the physical properties of z ^ 6 galax- 
ies. The key ingredient of these simulations is the inclu- 
sion of superwind feedback, critical to avoid the overpredic- 
tion by the simulations of the observed global star formation 
rate by reducing the reservoir of gas available for star forma- 
tion (Sprin gel & Hernquistll2003l) . Specifically, we used the 
"momentum-driven wind" model used inFinlatoret al. (2007, 
name ly, their "jvzw" model; see also lOppenheimer & Dave I 
|2006| for detailed descriptions). We have used the 32 h~' Mpc 
box simulation, combined with a 64 h~' Mpc box to better 
sample the bright end of the LF. We note that the key differ- 
ence between the AGN feedback implementation in the S AMs 
and the superwind feedback is that the former does not require 
star formation. 

Computing the SAM-predicted rest-frame V-band LFs is 
straightforward, as the catalog is complete in the luminosity 
range of interest and has no redshift errors^. We also extracted 
rest-frame colors of the galaxies in order to determine the LF 
for red and blue galaxies separately. 

4. RESULTS 

The comparison between the observed rest-frame V-band 
LFs of all galaxies at 2.7 < z < 3.3 and 2 < z < 2.5 with those 
predicted by the theoretical models is shown in Figure [T] It 
is immediately obvious that the models do not yet provide a 
precise description of galaxy evolution. Differences between 
the various models, and discrepancies between model predic- 
tions and data, are still as large as a factor of ~ 5 for certain 
luminosity and redshift ranges. 

At 2.7 < z < 3.3, the global LF predicted by the SAM of 
iBower et al. I fe)06) agrees well with the observed LF, al- 
though the SAM slightly underpredicts the density of galaxies 

- Available at http:// www.mpa-garching.mpg.de/Millennium] jLemsonI 
1200611 : see also footnote 6. 

^ To derive the model-predicted LF in a specific redshift interval, we av- 
eraged the number of galaxies as function of rest-frame V-band magnitude 
(with dust modeling included) of all redshift snapshots in the targeted red- 
shift interval. 
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Fig. 1 . — Compaiison between the rest-frame V-band observed global LFs 
and those predicted by models. The observed LFs are plotted with black cir- 
cles (1/V,i,ax method) with 1 a error bars (including lield-to-field variance) 
and by the black solid line (maximum likelihood method) with 1, 2, and 
3 (T solutions (gray shaded regions). The arrow shows the observed value 
of M*. Red lines show predictions f rom th ejBower et al. (2006) SAM, blue 
lines from the De Lucia & Blaizot] 120071) SAM, and green lines from the 
Finlator et al. (2007) SPH model. Poisson errors (1 a) are shown for the 
SPH model only, as they are very small for the SAMs. In the small panels, 
the ratio between the predicted and the observed LFs is plotted, together with 
the 1, 2, and 3 cr errors for the Bower et al. (2006) SAM (gray shaded re- 
gions). The oblique line regions delimit the comparison to the luminosity 
range probed by the sample o fjMarchesini et aL ilOOli) . 



around the knee of the LF. However, while at 2 < z < 2.5 the 
shape of the observed LF is broadly reproduced by the SAM, 
the predicted characteristic density <i>* is ~2.5 times larger 
than th e observed value. The SAM of iDe Lucia & Blaizot I 
(2007) has difficulty with both the normalization and the slope 
of the LF, which is too steep. At 2.7 < z < 3.3, the De Lucia 
model matches the faint end but underpredicts (by a factor of 
the bright end. At 2 < z < 2.5, instead, the predicted 
LF matches the bright end but overpred icts the faint end by a 
factor of > 2. The SPH simulations of Finlator et all (l2007h 
predict LFs that are qualitatively similar to those predicted 
by the two SAMs, although the former are characterized by 
larger uncertainties, due to the much smaller simulated vol- 
ume. 

We quantified these results by determining the luminosity 
density j\ (obtained by integrating the LF) for the observa- 
tions and models. The luminosity density is a more robust 
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Fig. 2. — Top panels: Ob.served luminosity density (j'^") as function of 
redshift of all (black circles), red (red circles), and blue (blue squares) galax- 
ies, splitting the sample based on rest-frame U — V (left panels) and B-V 
(right panels) colors. Bottom panels: Luminosity density predicted by the 
SAM of IBower et HMS) (j^^^) as function of redshift; symbols as in 
top panels; the observed evolution of j\ is also plotted with dashed lines for 
comparison. 



measure than M*, $*, and the faint-end slope a, because 
the errors in these parameters are highly correlated. The ob- 
served jy ( j^^) has been estimated by integrating the best-fit 
Schechter function down to My = -19.5, which is the faintest 
luminosity probed by the /T-selected sample'*. To estimate jy 
from the SAM (jy^'^), we have fitted the predicted LFs with a 
Schechter function, leaving M*, <!>*, and a as free parameters, 
applying the same limits as to the data. 

The comparison between and jy^ of iBower et al. I 
(120061) is shown in Figure|2](feoftom panels) by the black lines 
and data points. The Bower SAM matches the observed lumi- 
nosity density at z ^ 3. However, the model does not match 
the evolution of jy. In the model the luminosity increases 
with cosmic time, by a factor of ~ 1.6 from z ~ 3 to 2.2. 
By contrast, the observed luminosity density decreases with 
time, by a factor of ~ 1.8 over the same redshift range. Re- 
sults for the De Lucia SAM are similar, but for this model the 
difference between observed and predicted density is a strong 
function of the adopted faint-end integration limit. 

5. COLORS 

We investigated the cause of the discrepancies by splitting 
the sample into blue and red galaxies, using their rest-frame 
colors. Here we focus on the Bower model, as it provides the 
best match to the shape of the global LF, and a wide range 
of rest-frame colors are available. Interestingly, the results 
depend strongly on the choice of color: splitting the sample 
hy U -V color (as done in Marchesini et al. 2007) produces 
very different results than splitting hy B-V color 

To define red galaxies, we first use the criterion U - 

* As in 'Marchesini et al.' f2007) the 3 a error on jy was calculated by 
deriving the distribution of all the values of j\ within the 3 cr solutions of the 
Schechter LF parameters from the maximum-likelihood analysis, including in 
quadrature a 10% contribution from photometric redshift uncertainties. Using 
a brighter integration limit of the LF (My = -20.4) does not change the results 
of the comparison significantly. 



V >0.25, as done in iMarchesini et al.l (l2007b . As shown in 
the bottom left panel of Figure|2] the Bower model reproduces 
the densities of red and blue galaxies at z ~ 3 extremely well. 
The model overpredicts the densities of red and blue galaxies 
at z ~ 2.2, although it predicts the correct ratio between the 
two (roughly 1:1). 

Next, we use the criterion B-V >0.5.^ As can be seen in 
the top panels of Figured this criterion leads to very similar 
observed densities of red and blue galaxies as the U-V crite- 
rion. However, the predicted densities are in severe disagree- 
ment with the observations, particularly at z ~ 3 (see Fig. |2] 
bottom right panel). The red galaxy density at z ~ 3 under- 
predicts the observed density by a factor of 8. Qualita- 
tively similar results are obtained when jy^^ from the SAM 
of iDe Lucia & Blaizot I (l2007l) is used in the comparison.^ 

Irrespective of the color criterion that is used, we find that 
the predicted evolution of the red and blue luminosity densi- 
ties is in qualitative disagreement with the observed evolution. 
In the observations, the moderate evolution of the luminosity 
density is mainly driven by a decrease with cosmic time of the 
density of blu e galaxies, with the red galaxies evolving much 
less (see also IB rammer & van Dokkum Il2007h . By contrast, 
in the SAMs, the moderate evolution seen in the global LF is 
in the opposite sense and dominated by a strong evolution in 
the red galaxy population. 

6. DISCUSSION 

The main results of our comparison between the observed 
and the model-predicted res t-frame V-band LFs of galaxies 
at z > 2 are (1) the SAM of iBower et al. I (l2006l) reproduces 
well the observed LF at z ~ 3; (2) the models predict an in- 
crease with time of the rest-frame V-band luminosity density, 
whereas the observations show a decrease; (3) the models pre- 
dict strong evolution in the red galaxy population, whereas in 
the observations most of the evolution is in the blue popu- 
lation; (4) the models greatly underpredict the abundance of 
galaxies with B-V > 0.5 at z '--^ 3. 

The different results obtained for U-V and B-V colors 
are interesting, as they may hint at possible ways to improve 
the models. We further investigate the disagree ment between 
observ ed and predicted colors in the SAM of ' Bower et al. I 
( |2006|) in Figure [3] which shows the comparison of observa- 
tions and predictions in the B-V versus U-B diagram. While 
the SAM seems to broadly reproduce the observed U-B dis- 
tribution, it predicts galaxies that are systematically bluer in 
B-V than the observed galaxies. We have plotted evolution- 
ary trac ks of stellar population syn thesis models constructed 
with the'Bruzua l & Chariot I (l2003h code, assuming three dif- 
ferent prescriptions for the star formation history (SFH): a 
constant SFH (CSF), an exponentially declining in time SFH 
characterized by the parameter r = 300 Myr (t300), and an 
instantaneous burst model (SSP). We selected th e "Padova 
1994" evolutionary tracks, solar metallicity, the IChabrierl 
(I2003h initial mass function with lower and upper mass cut- 
offs 0.1 M,; and 100 M©, and modeled the e x tincti on by dust 
using the attenuation law of ICalzetti et al. I (12000 ). A new 
burst of star formation lasting 100 Myr and contributing 20% 
to the mass is also added at f = 2.1 x 10^ yr (f = 2.9 x 10^ yr) 
at z 3 (z ~ 2.2) to explore more complex SFHs. 

As can be deduced from Fig. [3] the differences between ob- 

^ For observed galaxies in the Marchesini et al. (2007) sample, U — V = 
0.25 roughly corresponds to B—V = 0.5. 

* The De Lucia model provides B — V colors, but no U-V colors. 
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Fig. 3. — B—VvsU-B comparison between observations {open circles) and predictions from the SAM of Bower et al. (2006, gray shaded regions) in tlie two 
targeted redshift intervals. Tlie error bars in tlie top left comer represent the median errors on the observed colors. The filled triangles with eiTor bars represent 
the mean colors of the observed galaxies and the error on the mean. The yellow, orange, and red lines show the evolutionary tracks described in §[6](CSF, t300, 
and SSP models, from top to bottom) with Av=0, 1, and 2, respectively. The tracks are plotted from 50 Myr to the age of the universe at the lower limit of the 
targeted redshift range. The cyan, blue, and purple lines show the evolution of the colors after a burst of star formation, for the three values of Ay, respectively. 
The arrow indicates the extinction vector for Av=l. The dashed lines correspond to B — 1^=0.5 and U-B = [U - V) - [B — V) = 0.25 - 0.5 = -0.25. Observed 
galaxies have redder B-V colors than predicted, possibly due to additional dust and/or secondary star bursts. 



served and predicted colors could be due to larger amount of 
dust and/or to more complex SFHs in the observed galaxies. 
The ad hoc treatment of dust absorption is a significant and 
well-known source of uncertainty in the models. Modifica- 
tions to the specific dust model could partly resolve the differ- 
ences between observations and SAM predictions. By simply 
multiplying the Ay in the SAM by a fixed factor, we were able 
to better reproduce the observed LFs at z ^ 2.2 (although mak- 
ing the faint-end slope of the red galaxy LFs quantitatively too 
steep) and to have a better agreement between observed and 
predicted colors. However, at z 3 this simple remedy is not 
able to solve the disagreement between the predicted and the 
observed number of B-V > 0.5 galaxies. We conclude that, 
while ad hoc modifications of the dust treatment might help to 
alleviate some of the found disagreements, it does not seem to 
be sufficient to accommodate the problem with global colors 
atz-3. 

While our ability to simulate galaxy formation has greatly 
improved in the past few years, our results imply that the 
present understanding of the physical processes at work in 



galaxy formation and evolution is still far from being satis- 
factory. On the observational side, more accurate redshift and 
color estimates would benefit studies of this kind. The lack of 
spectroscopic redshifts is particularly worrying, as systematic 
errors in redshift will lead to system atic errors in colors and 
luminosities (e.g.. lKriek et al. 1120061: M. Kriek et al. 2007, in 
preparation). 
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